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Abstract 

In early 2010, the Large Hadron Collider forward (LHCf) experiment mea- 
sured very forward neutral particle spectra in LHC proton-proton collisions. 
From a limited data set taken under the best beam conditions (low beam-gas 
background and low occurance of pile-up events), the single photon spectra at 
^/s=7^eV and pseudo-rapidity [i]) ranges from 8.81 to 8.99 and from 10.94 
to infinity were obtained for the first time and are reported in this paper. 
The spectra from two independent LHCf detectors are consistent with one 
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another and serve as a cross check of the data. The photon spectra are also 
compared with the predictions of several hadron interaction models that are 
used extensively for modeling ultra high energy cosmic ray showers. Despite 
conservative estimates for the systematic errors, none of the models agree 
perfectly with the measurements. A notable difference is found between the 
data and the DPMJET 3.04 and PYTHIA 8.145 hadron interaction models 
above 2 TeV where the models predict higher photon yield than the data. 
The QGSJET 11-03 model predicts overall lower photon yield than the data, 
especially above 2 TeV in the rapidity range 8.81<?7<8.99. 

Keywords: LHC, Ultra-High Energy Cosmic Ray, hadron interaction 
models 



1. Introduction 

The lack of knowledge about forward particle production in hadron colli- 
sions affects the interpretation of observations of Ultra-High Energy Cosmic- 
Rays (UHECR). Although UHECR observations have made notable improve- 
ments in the last few years [l| |2[ j3| 01 [sl 0] 01, some critical parts of the 
analysis depend on the Monte Carlo (MC) simulations of air shower devel- 
opment that are sensitive to the choice of the hadron interaction model. 
Accelerator data on the production of very forward emitted particles are in- 
dispensable for constraining the hadron interaction models but are usually 
not available from the large general purpose detectors. The Large Hadron 
Collider forward (LHCf) experiment has been designed to measure the neu- 
tral particle production cross sections at very forward collision angles of LHC 
proton-proton collisions, including zero degrees. When the LHC reaches its 
designed goal of 14 TeV collision energy, the energy in the equivalent labo- 
ratory frame will be 10^^ eV, a factor of one thousand increase compared to 
previous accelerator data in the very forward regions [sl 0]. 

Two detectors, called Arml and Arm2, have been installed in the in- 
strumentation slots of the TANs (Target Neutral Absorbers) located ±140 m 
from the ATLAS interaction point (IPl) and at zero degree collision angle. 
Inside a TAN the beam vacuum chamber makes a Y shaped transition from 
a single common beam tube facing the IP to two separate beam tubes joining 
to the arcs of LHC. Charged particles from the IP are swept aside by the inner 
beam separation dipole Dl before reaching the TAN so only neutral particles 
are incident on the LHCf detectors. This unique location covers the pseudo- 
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rapidity range from 8.7 (8.4 in case of the operation with the maximum beam 
crossing angle) to infinity (zero degrees). Each detector has two samphng 
and imaging calorimeters composed of 44 radiation lengths (1.55 hadron in- 
teraction lengths) of tungsten and 16 sampling layers of 3 mm thick plastic 
scintillators. The transverse sizes of the calorimeters are 20 mm x 20 mm and 
40mmx40mm in Arml, and 25mmx25mm and 32mmx32mm in Arm2. 
The smaller calorimeters cover the zero degree collision angle. Four X-Y 
layers of position sensitive detectors (scintillating fiber, SciFi, belts in Arml 
and silicon micro-strip sensors in Arm2; 1mm and 0.16 mm readout pitchs, 
respectively) are inserted in order to provide transverse positions of the show- 
ers. The LHCf detectors have energy and position resolutions for the elec- 
tromagnetic showers better than 5% and 200 fim, respectively, in the energy 
range >100GeV. More detail on the scientific goals, construction and per- 
formance of the detectors can be found in previous reports 11 1 12 1 13 



14 15 



This paper describes the first analysis results of LHCf data. Single pho- 
ton energy spectra are reported for i/i = 7 TeV proton-proton collisions. In 
SecE] the data set used in the analysis is introduced. In SecO] the analysis 
process and experimental results are presented. Beam related background 
and uncertainties are discussed in SecJH The experimental results are com- 
pared with MC predictions of several hadron interaction models in SeclDand 
summarized in SecEl 



2. Data 

Data used in this analysis was obtained on 15 May 2010 during proton- 
proton collisions at y/s=7TeV with zero degree beam crossing angle (LHC 
Fill 1104). The total luminosity of the three crossing bunches in this fill, 
L=(6. 3-6. 5) X 10^^ cm~^s~^, provided ideal operating conditions as discussed 
in SecJH The data that were taken during a luminosity optimization scan 
were eliminated from the analysis. The trigger for LHCf events was generated 
at three levels. The first level trigger (LIT) was generated from beam pickup 
signals (BPTX) when a bunch passed IPl. A shower trigger was generated 
when signals from any successive 3 scintillation layers in any calorimeter ex- 
ceeded a predefined threshold. Then the second level trigger for shower events 
(L2TA) was issued when the data acquisition system was armed. The thresh- 
old was chosen to achieve >99% efficiency for >100 GeV photons. Data were 
recorded with the third level trigger (L3T) when all the other types of second 
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level triggers (pedestal, laser calibration, etc) were combined. Examples of 
the longitudinal and lateral development of electromagnetic showers observed 
in the Arm2 detector are shown in FigJH In this case two electromagnetic 
showers from 7r° decay into two photons are shown, with each photon striking 
a different calorimeter of the Arm2 detector. The generation of the L2TA 
and L3T triggers, and hence the data recording, were performed indepen- 
dently for the Arml and Arm2 detectors. Data acquisition was carried out 
under 85.7% (Arml) and 67.0% (Arm2) average livetimes (eDAg)- The hve- 
times were defined as eoAQ = ^L2Ta/^ shower where N shower and Nl2ta are 
the number of counts in the shower and L2TA triggers, respectively. 

The integrated luminosities (J L dt) corresponding to the data used in 
this paper are 0.68 nb~^ (Arml) and 0.53 nb~^ (Arm2) after the data taking 
livetimes are taken into account. The absolute luminosity is derived from 
the counting rate of the Front Counters (FC) jTlj. FCs are thin plastic 
scintillators fixed in front of the LHCf main calorimeters and covering a wide 
aperture of 80mmx80mm. The calibration of the FC counting rates to the 
absolute luminosity was made during the Van der Meer scans on 26 April 
and 9 May 2010. The calibration factors obtained from the two scans differ 
by 2.1%. The estimated luminosities for the 15 May data differ by 2.7% 
between the two FCs. Considering the uncertainty of ±5.0% in the beam 



intensity measurement during the Van der Meer scans [16 
uncertainty of ±6.1% in the luminosity determination. 



we estimate an 



3. Analysis and Results 

3.1. Event Reconstruction 

The same analysis process has been adapted to each Arm independently. 
The transverse impact position of a particle is determined using the infor- 
mation provided by the position sensitive detectors. Using the position in- 
formation, the raw data from the scintillation layers are corrected for the 
non- uniformity of light collection and for shower 'leakage out' the edges of 



the layers [12] . Events that fall within 2 mm of the edges are removed from 
the analysis due to the large uncertainty in the energy determination owing 
to shower leakage. The recorded charge information is converted to deposited 
energy based on calibration runs with SPS fixed target experiments below 
200 GeV 10| 13j. The sum of the energy deposited in the 2"'^-13*'^ scintil- 



lation layers is converted to the primary photon energy by using a function 
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determined by MC simulation using the EPICS 8.81/COSMOS 7.49 simu- 



lation package [17[ and confirmed in the SPS beam tests. Note that this 
energy estimate does not represent the incident energy of hadrons because 
our calorimeters have only 1.55 hadron interaction lengths. In the detector 



simulations the Landau-Pomeranchuk-Migdal effect [18| [19| has been con- 
sidered and neglecting the LPM effect does not change the energy estimate 
at the 1% level because we sum the deposited energy up to a sufficiently deep 
layer. 

The linearity of each PMT was carefully tested before detector assembly 
over a wide range o f sig nal amplitude by exciting a scintillator using a 337 nm 
UV laser pulse l3| 1^ Although the measured non-linear response functions 



have been applied in the analysis, the difference between linear and non-linear 
reconstructions for 3TeV photons is only 0.5% at maximum. We also took 
data under LHC conditions with different PMT gains, but after applying gain 
calibrations no difference in the data sets was observed. Events having energy 
below 100 GeV are eliminated from the analysis to avoid corrections due to 
the trigger inefficiency and to reject particles produced in the interaction 
between collision products and the beam pipe. 

3.2. Single Event Selection 

To deduce the single photon energy, multi-hit events with more than one 
photons registered in a single calorimeter are eliminated. These multi-hit 
events are identified by using the lateral shower distribution measured by 
the position sensitive layers. According to MC simulation, the efficiency 
for correctly identifying true single photon events is >98%. The efficiency 
for identifying multi-hit events depend on the distance and the energy ratio 
of two photons and the detectors because of the different readout pitchs of 
the Arml SciFi belts and the Arm2 silicon micro-strip sensors. When the 
separation is greater than 1 mm and the lower energy photon has more than 
5% of the energy of the nearby photon, the efficiencies for identifying multi- 
hit events are >70% and >90% for Arml and Arm2, respectively. 

To estimate the systematic uncertainty in the multi-hit identification ef- 
ficiency, we produced an artificial sample of multi-hit event sets by superim- 
posing two clearly single photon-like events for both the experimental and 
MC data based on the EPOS 1.99 model. Details of the MC simulations are 
described in SecEl To choose the energies and separation of a photon pair, 
we followed the distributions determined by the DPMJET 3.04 model. For 
the two artificial data sets the efficiencies for identifying multi-hit events do 
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not differ by more than 10% and 3% over the entire energy range for Arml 
and Arm2, respectively. This affects the final single photon energy spectrum 
shape by less than 1% below 1.5 TeV and increasingly up to 2-20% at 3 TeV. 
Tha maximum difference is found for the Arml large calorimeter. 

Next, we compared the effect of the multi-hit cut on the Arml and Arm2 
detectors. While the fraction of events thrown out by this cut differs by less 
than 5% between 0.5 TeV and 1.5 TeV, it gradually increases to 30% and 
60% at 3 TeV for the small and large calorimeters, respectively. The main 
reason for these differences is the different geometry of the Arml and Arm2 
calorimeters. The different performances of the position sensitive detectors 
in the two Arms and an uncertainty in the absolute energy scale discussed 
in Sec J3.4l may also contribute to the differences in multi-hit identification 
fractions of the Arml and Arm2 detectors. Because we cannot presently 
separate the sources of the difference and hence cannot apply corrections to 
the data, we assign the differences divided by \/2 as part of the systematic 
uncertainty for each detector. Finally we take quadratic sum of the two 
uncertainties related to the multi-hit identification efficiency and the multi- 
hit cut as systematic error of the single photon selection procedure. 

3.3. Photon Event Selection 

To select only electromagnetic showers and eliminate hadron (predomi- 
nantly neutron) contamination, a simple parameter, L9o% is defined. L9o% is 
the longitudinal distance in radiation lengths measured from the entrance to 
a calorimeter to the position where 90% of the total shower energy has been 
deposited. Figj2] shows the distribution of L9o% for the 20 mm calorimeter 
of the Arml detector for the events with the reconstructed energy between 
500 GeV and 1 TeV. Two distinct peaks are observed corresponding to photon 
and hadron (neutron) events. The L9o% distributions for pure photon and 
hadron samples are generated by MC simulation using the collision prod- 
uct generator QGSJET 11-03 as shown in Figj2l They are called 'templates' 
hereafter. The choice of hadron interaction model in determining the tem- 
plate does not affect the results in this paper. In the event selection, we set 
an energy dependent criteria in L9o% to keep the photon detection efficiency 
ep/D=90% over the entire energy range based on the photon template. The 
purity of the selected photon events is determined by normalizing the tem- 
plate functions to the observed L9o% distribution. The purity, P, is defined 
as P=Npftof/ {^phot+^had) in each energy bin. Here ^phot and are the 
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numbers of photon and hadron (neutron) events in the templates in the se- 
lected L9o% range. Multiplying each energy bin by PxeJ,]^, we obtained non 
biased photon energy spectra. 

Some disagreements in the Lqq% distribution are found between the data 
and the MC calculations. This may be caused by errors in the absolute 
energy determination and channel-to-channel calibrations and may also be 
motivation for studying the LPM effect in detail. Here we consider a system- 
atic uncertainty caused by the uncertainty of the template fitting method in 
the correction of the photon spectra. Small modifications of the template 
functions, widening with respect to the peak position up to 20% and con- 
stant shift up to 0.7 radiation lengths, to give the best match with the data, 
provide another estimate of the correction to the photon spectra. The differ- 
ence of the correction factors between the original and the modified template 
methods amount to 5-20% from low to high photon energy and this is as- 
signed as a systematic uncertainty of the particle identification in the final 
spectra. 

3.4- Energy Scale Uncertainty from Mass Reconstruction 

When each of two calorimeters records a single photon as shown in FigJH 
shower 'leakage in' is corrected according to a function based on MC sim- 
ulation. Using the corrected energies and positions of the shower axes, the 
invariant mass of the photon pair is calculated assuming their vertex is at the 
interaction point. In MC simulations of the full detector response and the 
analysis process, we confirmed the reconstructed mass peaks at 135.2 MeV 
in Arml and 135.0 MeV in Arm2, thus reproducing the 7r° mass. The statis- 
tical uncertainty in the reconstructed invariant mass of the MC simulations 
is ±0.2 MeV. 

On the other hand, the reconstructed invariant masses of photon pairs 
for the experimental data are 145.8±0.1MeV (Arml) and 140.0±0.1MeV 
(Arm2) where ±0.1 MeV uncertainties are statistical. A portion of the 7.8% 
and 3.7% invariant mass excess compared to the 7r° mass reconstructed in 
the MC simulations can be explained by the well understood systematic error 
of the absolute energy scale, estimated to be ±3.5%. This 3.5% systematic 
error is dominated by the errors in factors converting measured charge to 
deposited energy and by the errors in corrections for non-uniform light col- 
lection efficiency. Uncertainties in determining the opening angle of a photon 
pair and the shower leakage-in correction, typically ±1% and ±2% respec- 
tively, are also sources of error in mass reconstruction. These known elements 
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quadratically add up to a systematic mass shift of 4.2% and can explain the 
mass shift in the Arm2 detector, but not Arml. 

Because all the two photon invariant mass shift may not be due to the 
energy scale uncertainty, we did not apply any correction for energy scale 
in the energy spectra presented in Sec J3.5[ Instead we assigned asymmetric 
systematic errors in the absolute energy scale, [-9.8%, +1.8%] (Arml) and 
[-6.6%, +2.2%] (Arm2). Here we assumed uniform and Gaussian probability 
distributions for the energy scale errors estimated from the mass shift (7.8% 
and 3.7%) and the known systematics (3.5%), respectively. After the stan- 
dard deviations of two components (7.8%/-\/3 in case of the uniform prob- 
ability distribution) were quadratically added, the systematic error bands 
are assigned with respect to the central value of the mass shift. To deter- 
mine the systematic errors in the final energy spectra, we reconstructed two 
energy spectra by scaling the energy using the two extremes quoted above. 
The differences from the non-scaled spectrum to the two extreme spectra are 
assigned as systematic errors in each energy bin. 

3.5. Spectra Reconstruction 

To compensate for the different geometry of two arms, we selected com- 
mon rapidity and azimuthal ranges to deduce the photon energy spectra. 
The ranges used for the small calorimeters and the large calorimeters are 
[r;>10.94, A(^=360.0°] and [8.99>r]>8.81, A0=2O.O°], respectively. Here r/, 
and A(/) represent the pseudo-rapidity, azimuthal direction and interval of 
with respect to the beam axis which is centered on the small calorimeters. 

Photon spectra measured in the small and large calorimeters are shown 
in Figj3l The red and blue plots show the results from Arml and Arm2, 
respectively. The error bars and shaded areas indicate the one standard 
deviation statistical and systematic errors, respectively, uncorrelated between 
the two detectors. On the vertical axis, the number of inelastic collisions, 
Nj„e, is calculated as Nj„e = c"me /Ldt assuming the inelastic cross section 
o"me = 71.5 mb. Using the integrated luminosities introduced in Secj2l Nj„e 
= 4.9x10^ for Arml and 3.8x10'' for Arm2. From Fig 0] we find general 
agreement between the two Arms that are within the errors. The reason 
for the difference between the two Arms in the small calorimeters (higher 
rapidity) is not yet understood. However, because the difference is still within 
the errors, we did not apply any correction. 

Here we note that the obtained spectra are expected to be distorted from 
the single photon inclusive spectra due to the analysis processes especially 
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the multi-hit cut that reflects the physical size of the LHCf detector and 
the difl^erences between the hadron interaction models. The multi-hit cut is 
expected to suppress the event rate per pp interaction while the inefficiency 
of multi-hit identification raises the event rate at high energy because we 
misreconstruct multi photon energy as single photon energy. To estimate 
the deformation of the energy spectra the reconstructed energy spectra of 
photons normalized to the true inclusive single photon spectra are studied 
by MC simulations for the different hadron interaction models. Here in cal- 
culating the 'true inclusive single photon spectra' particle decay in the 140 m 
flight path from the interaction point to the LHCf detectors and the LHCf 
calorimeter aperture are taken into account. In the case of multi- hits in a sin- 
gle calorimeter, each photon is counted independently. As a result, 0-15% of 
energy independent supression is found below 2 TeV and it turns to gradually 
rise up to +15% over 3 TeV. The maximum difference between interaction 
models and between the two arms are about 10% and 5%, respectively. 

4. Beam Related Background and Uncertainties 

The events containing more than one collision (pile-up events) in a single 
bunch crossing may cause an additional bias. Given that a collision has 
occurred, the probability of pile-up (P(n>2)/P(n>l)) can be calculated from 
the Poisson probability distribution. Using the highest bunch luminosity 
of L=2.3xl0^^ cm~^s~^ used in this analysis, inelastic cross section (7j„e = 
71.5 mb and the revolution frequency of LHC frev — 11-2 kHz, the probability 
is P(n>2) /P(n>l)=0.072. Considering the acceptance of a LHCf calorimeter 
for an inelastic collision, ~0.03, only 0.2% of events have more than one event 
due to the pile- up and they are eliminated in the multi- hit cut. We conclude 
that pile-up does not affect our analysis. 

In the geometrical analysis of the data, we assumed the projected position 
of the zero degree coUision angle at the LHCf detectors, referred to as the 
'beam-center' hereafter, can move from flU to fill owing to slightly different 
beam transverse position and crossing angles at the IP. We determined the 
'beam center' at the LHCf detectors by two methods; first by using the 
distribution of particle impact positions measured by the LHCf detectors and 
second by using the information from the Beam Position Monitors (BPMSW) 
installed ±21 m from the IP. From the analysis of the fills 1089-1134, we 
found a maximum ~4mm shift of the beam center at the LHCf detectors, 
corresponding to a crossing angle of ~30 /xrad assuming the beam transverse 
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position did not change. The two analyses gave consistent results for the 
location of the beam center on the detectors within 1 mm accuracy. In the 
geometrical construction of events we used the beam-center determined by 
LHCf data. We derived photon energy spectra by shifting the beam-center 
by 1 mm. The spectra are modified by 5-20% depending on the energy and 
the rapidity range. This is assigned as a part of systematic uncertainty in 
the final energy spectra. 

The background from collisions between the beam and the residual gas in 
the vacuum beam pipe can be estimated from the data. During LHC oper- 
ation, there were always bunches that did not have a colliding bunch in the 
opposite beam at IPl. We call these bunches 'non-crossing bunches' while the 
normal bunches are called as 'crossing bunches.' The events associated with 
the non-crossing bunches are purely from the beam-gas background while the 
events with the crossing bunches are mixture of beam-beam collisions and 
beam-gas background. Because the event rate of the beam-gas background 
is proportional to the bunch intensity, we can calculate the background spec- 
trum contained in the crossing bunch data by scaling the non-crossing bunch 
events. We found the contamination from the beam-gas background in the 
final energy spectrum is only ~0.1%. In addition the shape of the energy 
spectrum of beam-gas events is similar to that of beam-beam events, so 
beam-gas events do not have any significant impact on the beam-beam event 
spectrum 

The collision products and beam halo particles can hit the beam pipe and 
produce particles that enter the LHCf detectors. However according to MC 



simulations, these particles have energy below 100 GeV [10[ and do not affect 
the analysis presented in this paper. 

5. Comparison with Models 

In the top panels of FigJH photon spectra predicted by MC simulations 
using different modeK QGSJET 11-03 (blue) jli, DPMJET 3.04 (red) |2l|, 
SIBYLL 2.1 (green) EPOS 1.99j(magenta)|23 and PYTHIA 8.145 (de- 



fault parameter set; yellow) [2J] [25[ for collisions products are presented 



together with the combined experimental results. To combine the experi- 
mental data of the Arml and Arm2 detectors, the content in each energy 
bin was averaged with weights by the inverse of errors. The systematic un- 
certainties due to the multi-hit cut, particle identification (PID), absolute 
energy scale and beam center uncertainty are quadratically added in each 
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energy bin and shown as gray shaded areas in FigjH The uncertainty in the 
luminosity determination (±6.1% as discussed in SecJ2]), that is not shown 
in FigJU can make an energy independent shift of all spectra. 

In the MC simulations, 1.0x10^ inelastic collisions were generated and 
the secondary particles transported in the beam pipe. Deflection of charged 
particles by the Dl beam separation dipole, particle decay and particle in- 
teraction with the beam pipe are taken into account. The responses of the 
detectors were calculated using the EPICS/COSMOS libraries taking the 
random fluctuation equivalent to electrical noise into account. The same 
analysis procedures were then applied to the MC simulations as to the ex- 
perimental data except for the particle identification and its correction. In 
the analysis of the MC data set, we used the known information of the particle 
type. In the bottom panels the ratios of MC simulations to the experimental 
data are plotted together with the statistical and systematic uncertainties of 
the experimental data. The statistical uncertainty of the EPOS 1.99 model 
is also plotted as magenta shaded representative of the various 

models. 

We find that none of the models lies within the errors of our data over 
the entire energy range. Some remarkable features are: 

1) DPMJET 3.04 and PYTHIA 8.145 show very good agreement with 
the experimental result between 0.5 and 1.5 TeV for ?7>10.94, but they pre- 
dict significantly larger photon yield at high energy >2 TeV in both rapidity 
ranges. 

2) QGSJET 11-03 predicts overall lower photon yield than the experimen- 
tal result. This is significant above 2 TeV in the rapidity range 8.81<r7<8.99. 

3) For 7]>10.94, SIBYLL 2.1 shows a very good agreement with the ex- 
perimental result for the spectral shape for >0.5TeV, but predicts a photon 
yield only half of the experimental result over the entire energy range. 

6. Summary 

LHCf has measured for the first time the single photon energy spectra of 
high energy photons in the very forward region of proton-proton collisions at 
LHC. After selecting data with common rapidity ranges, the two indepen- 
dent LHCf detectors (Arml and Arm2) installed on either side of IPl gave 
consistent results in two rapidity ranges even though the geometrical accep- 
tances of the two detectors differ. The combined spectra of the two detectors 
are compared with the prediction of various hadron interaction models. It 
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is found that none of the model predictions have perfect agreement with the 
experimental results within statistical and systematic errors. 

The conservative systematic errors assigned in this analysis will be im- 
proved upon in future studies. In addition studies of other measurements 
like the inclusive single photon spectra, inclusive single vr" and neutron pro- 
duction spectra, and neutral particle transverse momentum spectra are now 
ongoing using the already accumulated LHCf data. By combining the LHCf 
data with the recent studies on particle production in the central rapidity 



region of LHC collisions [26j it is now for the first time possible to make 
critical tests of hadron interaction models by using collider data over a very 
wide rapidity range. 
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Figure 1: An example of the tt*^ candidate events observed in the LHCf Arni2 detector. 
Top two figures show the longitudinal developments of the two photon initiated showers 
observed in the 25 mm and the 32 mm calorimeters. In the middle and bottom panels, 
transverse X and Y profiles of the showers are shown. Different colors indicate data in the 
different silicon layers. 
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Figure 2: The Lgg% distribution measured by the Arml 20 mm calorimeter for the re- 
constructed energy of 500 GeV - 1 TeV. Plots are experimental data and the red and 
blue histograms are the templates calculated from the pure photon and pure hadron MC 
events, respectively. The two templates are independently normalized to best describe the 
observed data. 
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Figure 3: Single photon spectra measured by the Armi (red) and Arm2 (blue) detectors. 
Left (right) panel shows the results for the small (large) calorimeter or large (small) rapid- 
ity range. The error bars and shaded areas indicate the statistical and systematic errors, 
respectively. To discuss consistency of two detectors, only uncorrelated components are 
plotted for the systematic errors. 
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Figure 4: Comparison of the single photon energy spectra between the experimental data 
and the MC predictions. Top panels show the spectra and the bottom panels show the 
ratios of MC results to experimental data. Left (right) panel shows the results for the large 
(small) rapidity range. Different colors show the results from experimental data (black), 
QGSJET 11-03 (blue), DPMJET 3.04 (red), SIBYLL 2.1 (green), EPOS 1.99 (magenta) 
and PYTHIA 8.145 (yellow). Error bars and gray shaded areas in each plot indicate the 
experimental statistical and the systematic errors, respectively. The magenta shaded area 
indicates the statistical error of the MC data set using EPOS 1.99 as a representative of 
the other models. 
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